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Executive  Summary 


Ship  design  practice  has  been  to  measure  stability  by  static  criteria  and  to  compensate  for  dynamic 
effects  through  a  margin  of  safety.  However,  there  is  a  fundamental  difference  between  static  and 
dynamic  stability.  The  Fastnet  Yacht  Race  disaster  in  1979  revealed  that  "certain  factors  which 
result  in  favorable  static  stability  characteristics  may  actually  present  greater  danger  when 
considered  in  light  of  a  dynamic  analysis"  (Stephens,  Kirkman  and  Peterson,  1981). 

The  existing  linear  strip-theory  method  cannot  be  used  for  assessing  capsizing.  Advanced  nonlinear 
simulation  methods  are  required.  As  we  shall  see,  such  advanced  methods  are  now  under 
development  and  their  application  to  the  assessment  of  the  vessel  dynamic  stability  problem  is  a 
realistic  practical  goal  today. 

The  main  objective  of  the  present  project  has  been  to  investigate  the  capabilities  of  the  3-D 
nonlinear  time-domain  Large- Amplitude  Motion  Program  (LAMP)  for  the  evaluation  of  fishing 
vessels  operating  in  extreme  waves.  The  project’s  focus  was  building  upon  the  previous  LAMP 
development  and  extending  it  to  the  modeling  of  maritime  casualties,  including  a  time-domain 
simulation  of  a  ship  capsizing  in  stem  quartering  seas.  This  modeling  capability  will  allow  both  the 
analysis  of  recorded  casualties  and  the  identification  of  potential  safety  concerns. 

Ship  motions  in  stem  quartering  seas  are  extremely  complicated  since  the  roll  motion  is  highly 
nonlinear  and  the  viscous  effects  are  important.  A  typical  example  of  stem  quartering  sea  capsizing 
is  illustrated  in  a  time  sequence  in  Figure  1 .  The  simulations  are  for  the  70-foot  (2 1 .3  m)  fishing 
vessel  Italian  Gold  with  a  heavy  loading  condition  in  regular  stem  quartering  seas.  The  center  of 
gravity,  CG,  of  the  ship  is  located  at  the  mid  ship  and  0.27  ft  (0.08  m)  above  the  load  waterline.  A 
linear  moderate  wave  (wave  height,  h=6  ft  [1.8  m],  and  wavelength  X=100  ft  [30.5  m])  approaches 
from  the  stem  quarter  in  the  starboard  direction.  The  simulation  shows  that  operating  with  a  course¬ 
keeping  autopilot,  the  ship  turns  into  almost  beam  sea  condition  and  capsizes  in  the  direction  of 
wave  propagation. 

In  this  report,  the  methodology  used  to  assess  vessel  stability  and  safety  is  discussed.  It  emphasizes 
the  importance  of  performing  analyses  of  dynamic  stability  rather  than  applying  margins  of  safety 
to  static  stability  criteria.  The  need  for  computational  tools  for  accurate  stability  and  safety 
assessment  is  addressed. 

Using  the  current  LAMP  code,  the  present  study  shows  some  examples  of  fishing  boat  capsizings 
in  seas  that  are  initially  off  the  stem  quarter.  These  results  clearly  demonstrate  the  necessity  and 
power  of  a  nonlinear  time-domain  simulation  tool  for  the  study  of  vessel  stability  and  for  the 
assessment  of  ship  safety.  However,  an  extensive  validation  study  and  possibly  further 
improvements  to  the  present  method  may  be  required  for  accurate  predictions  of  extreme  ship 
motions  in  more  general  extreme  sea  conditions. 
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1  Introduction 


One  of  the  primary  missions  of  the  U.S.  Coast  CJuard  is  the  protection  of  life  and  property  by  the 
establishment  and  enforcement  of  marine  safety  standards.  The  ideal  standard  seeks  to  ensure 
safety  without  unduly  affecting  the  ship's  operability.  A  major  safety  concern  is  the  prevention  of 
the  loss  of  life  and  property  due  to  ship  capsizing.  Current  regulations  seek  to  prevent  such 
occurrences  by  setting  minimum  stability  and  freeboard  requirements.  These  regulations  are  based 
mainly  on  hydrostatics.  They  were  developed  from  an  analysis  of,  and  experience  with,  traditional 
ship  configurations. 

The  stability  assessment  of  new,  innovative  ship  forms  and  the  assessment  of  capsizing  accidents 
often  require  very  expensive  and  time  consuming  experiments.  The  existing  ship  motion  prediction 
tools  are  primarily  based  on  hydrostatics  and  linear  strip  theory,  which  can  only  be  used  for 
assessing  small  amplitude  motions  in  moderate  sea  conditions.  Therefore,  an  accurate  ship  motion 
simulation  method  may  have  a  large  impact  on  ship  safety  assessment. 

Computational  simulation  techniques  and  computer  architectures  have  finally  reached  such  a  level 
of  sophistication  that  the  development  of  a  simulation  system  for  vessel  stability  and  safety 
assessment  for  extreme  seas  is  a  practical  goal.  The  purpose  of  this  report  is  to  discuss  the  recent 
advances  in  computational  hydrodynamics  research  and  the  related  practical  engineering  systems, 
in  particular  the  LAMP  System  (Lin  and  Yue,  1990,  Lin,  et  al,  1992,  1993,  1994),  for  the 
assessment  of  the  stability  and  safety  of  a  vessel  operating  in  extreme  seas. 

Since  stability  criteria  are  primarily  based  on  static  stability,  it  is  extremely  important  to  emphasize 
that  the  physics  governing  static  stability  is  quite  different  from  the  physics  for  dynamic  stability. 
To  illustrate  this  point,  we  shall  first  look  at  a  sailing  yacht  disaster  which  has  been  investigated 
extensively  and  which  is  quite  well  understood. 


1.1  Dynamic  vs.  Staf/c  Stability 

The  Fastnet  Race  of  1979  is  considered  to  be  the  greatest  disaster  in  the  history  of  the  sport  of 
yachting.  Seventy-seven  boats  were  completely  capsized  and  fifteen  sailors  died  (Rousmaniere, 
1980).  Stephens,  Kirkman  and  Peterson  (1981)  analyzed  the  Fastnet  disaster  in  their  landmark 
paper  on  "Sailing  Yacht  Capsizing".  They  addressed  the  capsize  mechanism,  the  environmental 
conditions,  and  the  design  approach  which  led  to  the  terrible  disaster.  They  pointed  out  that  design 
practice  has  been  to  measure  stability  by  static  criteria  and  to  compensate  for  dynamic  effects 
through  margins  of  safety.  Their  investigation  of  the  Fastnet  Race  disaster  revealed  that "  certain 
factors  which  result  in  favorable  static  stability  characteristics  may  actually  present  greater  danger 
when  considered  in  light  of  a  dynamic  analysis." 

This  is  a  very  important  aspect  of  vessel  stability  that  unfortunately  is  often  overlooked  in  setting 
safety  requirements.  For  example,  it  is  often  assumed  that  a  vessel's  stability  is  a  function  of  its 
freeboard  with  larger  freeboard  providing  greater  capsizing  resistance.  This  is  correct  from  a  static 
point  of  view,  but  Stephens,  etal,  showed  that  the  dynamics  of  the  single  wave  impact  capsizing 
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mechanism  which  dominated  the  Fastnet  79  casualties  (see  Figure  2)  had  the  opposite  effect.  They 
asserted  that; 

. .  .freeboard,  which  helps  raise  the  zero-stability  crossing  in  a  static  case 
and  hence  appears  as  safe,  is  the  source  of  much  overturning  energy  being 
impacted  to  the  yacht  due  to  the  large  area  being  struck  by  the  breaker  and 
the  increased  moment  arm  acting  for  overturn. 


Stage  1  - 
Approach  of 
Breaking  Wave 


Stage  2  - 
Hull  Response 
to  Wave  Slope 


Stage  3  - 
Breaking 
Wave  Impact 


Stage  4  - 
Capsize 


Figure  2.  Schematic  of  the  Single  Wave  Impact  Capsizing  that  Dominated  the 
Fastnet  ’79  Casualties  (from  Stephens,  et  al,  1981) 


Furthermore,  Stephens,  et  al.,  concluded  that 

. .  .the  beam  contribution  to  static  stability  is  washed  out  in  a  capsize  by  a 

corresponding  moment  caused  by  local  wave  slope. 

The  reason  for  reviewing  the  Fastnet  Race  disaster  is  to  stress  the  importance  of  analyzing  dynamic 
stability  in  extreme  sea  conditions  and  to  focus  attention  on  the  fact  that  there  is  a  fundamental 
difference  between  static  and  dynamic  stability.  Vessel  safety  requirements  cannot  be  established 
by  considering  static  stability  alone  and  then  applying  some  safety  factor  to  include  dynamic 
aspects.  Design  trends  driven  by  static  stability  requirements  resulted  in  a  catastrophic  disaster  in 
the  Fastnet  Race  of  1979  because  they  increased  freeboard  and  increased  beam,  causing  dynamic 
stability  problems. 
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1.2  Summary  of  Present  Study 

In  this  report,  we  address  the  research,  development,  and  application  of  advanced  computational 
methods  for  the  assessment  of  dynamic  stability  in  a  seaway.  As  an  example,  we  use  the  LAMP 
system  to  study  a  particular  casualty:  the  capsize  of  the  stem  dragger  Italian  Gold  in  a  storm  off 
Massachusetts  in  September  of  1994.  The  intent  of  this  exercise  was  two-fold  : 

•  To  attempt  to  leam  more  about  the  dynamic  mechanisms  causing  capsizing  in  statically 
stable  vessels. 

•  To  test  and  improve  the  LAMP  system  for  application  to  such  casualty  analysis. 

Two  features  which  were  recently  added  to  LAMP  have  made  it  possible  to  do  these  analyses  in 
other  than  head-seas  conditions:  the  ability  to  calculate  viscous  forces  including  roll  and  appendage 
lift,  and  a  dynamic  automatic  control  system  for  steering.  This  challenging  real-world  problem 
provided  an  opportunity  to  test  these  features  and  improve  them. 

Static  stability  is  assumed  to  be  well  imderstood.  It  is  important  to  recognize  that  the  advanced 
computational  hydrodynamic  tools  discussed  in  this  report  are  equally  important  for  addressing  all 
of  the  following  problems  related  to  vessel  safety  in  waves: 

•  Dynamic  Stability  Stractural 

•  Equipment  Damage 

•  Crew  Safety. 

However,  in  this  report,  we  focus  on  the  dynamic  stability  problem.  Stability  as  it  is  discussed  here 
will  include  both  intact  and  damaged  stability. 

The  development  of  computational  tools  for  dynamic  vessel  stability  must  be  considered  as  a 
portion  of  the  much  larger  topic  of  seakeeping.  Seakeeping  assessment,  which  includes  both  the 
wave-induced  motions  and  hydrodynamic  loads,  may  be  divided  into  two  classes: 

•  Linear  frequency-domain  predictions 

•  Nonlinear  time-domain  simulations. 

Linear  frequency-domain  prediction  methods  {i.e.  Strip  Theory)  have  been  extremely  successful  in 
determining  sea  state  operability  limitations  for  weapon  systems  on  naval  vessels  (Kennel,  1985). 
Such  methods  have  also  been  useful  in  estimating  the  wave-induced  loads  for  large  ships  (Liu,  et 
ai,  1992).  However,  the  linear  "strip  theory"  tools  are  based  upon  the  assumption  that  both  the 
motions  and  the  wave  amplitudes  are  small  relative  to  the  vessel’s  dimensions  (in  particular  the 
draft).  Furthermore,  many  linear  methods  assume  wall  sidedness.  These  are  serious  limitations  -  the 
assumptions  are  not  valid  in  general  for  vessel  response  in  extreme  seas. 

Nonlinear  time-domain  simulation  is  required  to  determine  the  vessel’s  response  in  extreme  seas. 
Because  of  this  requirement,  dynamic  stability  predictions  are  an  order  of  magnitude  more  complex 
than  linear  frequency-domain  predictions.  For  example,  the  wave  field  description  for  extreme 
response  prediction  must  contain  much  more  detailed  information  than  the  wave  energy-spectrum 
representation  used  for  linear  prediction. 
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In  order  to  obtain  the  probabilistic  estimates  needed  for  setting  safety  standards,  one  has  to  apply  a 
combination  of  deterministic  and  probabilistic  calculations.  The  assessment  of  a  vessel’s  dynamic 
stability  in  waves  may  be  divided  into  three  parts: 

•  Wave-Event  Modeling:  extreme  wave  characterization,  selection  of  potentially  dangerous 
extreme  wave  events,  and  detailed  numerical  modeling  of  the  complex  nonlinear 
hydrodynamics  aspects  of  the  selected  wave  events 

•  Vessel-Response  Simulations:  an  accurate  time-domain  simulation  of  the  vessel’s 
response  to  the  selected  wave  events 

•  Probabilistic  Predictions:  an  estimate  of  the  probability  of  occurrence  of  the  wave/vessel 
encounters  that  will  result  in  catastrophic  responses. 

The  importance  and  development  of  these  three  parts  have  been  addressed  by  Salvesen  and  Lin 
(1993)  in  their  proposed  SAFE  SEAS  System.  TTie  first  and  the  third  parts  will  not  be  discussed 
further  in  this  report.  The  current  development  of  the  vessel-response  simulations  will  be  discussed 
here. 

Section  2  of  this  report  gives  a  general  description  of  the  LAMP  System  used  for  the  current 
simulation  study.  Section  3  describes  the  hull  form  studied  and  the  hydrostatic  characteristics  of  the 
vessel.  Results  of  using  the  LAMP  system  in  nonlinear  simulation  of  ship  capsizing  in  stem 
quartering  seas  are  presented  in  Section  4.  Both  time-domain  simulations  and  some  of  the 
mechanism  that  causes  the  ship  to  capsize  in  stem  quartering  seas  will  be  given. 
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2  The  LAMP  System 

2.1  The  LAMP  Approach 

In  1990,  Lin  and  Yue  presented  a  three-dimensional  time-domain  method  to  study  large-amplitude 
motions  and  loads  of  floating  bodies  in  waves.  In  their  so-called  "body-exact"  approach,  the  free- 
surface  boundary  conditions  are  linearized  and  the  body  boundary  condition  is  satisfied  exactly  on 
the  portion  of  the  instantaneous  wetted  surface  that  lies  below  the  undisturbed  free  surface.  The 
problem  is  solved  using  a  transient  free-surface  Grreen's  function  singularity  distribution.  The 
validity  and  practical  utility  of  this  method  has  been  demonstrated  by  several  studies,  including 
predictions  of  large-amplitude  motion  coefficients,  motion  history  of  a  ship  advancing  in  an 
irregular  seaway,  and  the  effect  of  bow  flare  on  wave  loads  (see  Lin  and  Yue  1990,  1992;  Lin  et  al., 
1991,  1992). 

In  1993,  Lin  and  Yue  extend  the  applicability  of  their  method  to  allow  ship  motions  in  more  severe 
wave  conditions  in  which  both  the  body  motions  and  the  incident  waves  can  be  large.  In  this  new 
Large-Amplitude  Motion  Program,  LAMP,  the  body  boimdary  condition  is  satisfied  on  the 
instantaneous  wetted  surface  below  the  incident  wave  profile  with  the  assumption  that  the 
diffracted  waves  are  small  compared  to  the  incident  wave  and  that  the  incident  wave  slopes  are 
small.  At  each  time  step,  local  incident  free  surface  elevations  are  used  to  transform  the  body 
geometry  into  a  computational  domain  with  a  deformed  body  and  a  flat  free  surface.  By  linearizing 
the  free  surface  boundary  conditions  about  this  incident  wave  surface,  the  problem  can  be  solved  in 
the  computational  domain  using  linearized  free-surface  transient  Green's  fimctions.  The  two  main 
features  of  this  new  large-amplitude  approach  are: 

1 .  true  hydrodynamic  effects  for  the  wetted  portion  of  the  ship  under  the  incident  wave  surface 

2.  automatic  inclusion  of  the  correct  hydrostatic  and  Froude-Krylov  forces. 

A  summary  of  this  new  approach  is  given  in  Lin  et  al.  (1994). 

The  LAMP  approach  solves  the  six-degree-of-freedom  d5mamical  equations  of  motion  in  the  time 
domain  for  motion  simulation.  At  each  time  step,  forces  and  moments  acting  on  the  body  are  given. 
For  solving  ship  motion  problems  in  oblique  seas,  viscous  and  lifring  forces  are  important.  A  brief 
description  of  the  viscous  force  calculations  included  in  the  LAMP  approach  is  given  in  the  next 
section. 

In  the  LAMP  code,  the  seaway  is  represented  by  summation  of  linear  wave  components.  The 
seaway  can  be  specified  as  either  regular  (single  component)  or  a  random  wave  (multiple 
components)  specified  by  a  wave  spectrum.  A  discussion  of  the  seaway  representation  is  given  later 
in  this  section.  For  motion  simulation  of  fishing  vessel  in  stem  quartering  seas,  course-keeping 
mdder  control  is  used.  A  simple  proportional,  integral,  and  derivative  (PID)  control  of  the  mdder  is 
implemented  in  the  LAMP  code.  This  control  algorithm  is  discussed  briefly  also  later  in  the 
section. 
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2.2  Viscous  Forces  in  the  Time  Domain 


2.2.1  Overview 

In  oblique  or  beam  seas,  forces  due  to  viscous  and  lift  effects  will  have  a  significant  effect  on  the 
motions  and  loads.  LAMP  includes  an  option  to  approximate  some  of  these  effects  in  the  time- 
domain.  These  included  effects  may  not  completely  characterize  the  viscous  flow  separation  effects 
present  in  maneuvering  cases.  Further  development  using  full  viscous  flow  codes  such  as  Re5molds 
Averaged  Navier  Stokes  (RANS)  solvers  may  be  necessary  to  better  characterize  the  3-D  flow 
effects  found  in  maneuvering  cases.  The  viscous  and  lift  effects  approximated  are  as  shown  in 
Table  1 .  For  each  effect,  the  table  presents  a  reference  for  the  calculation  method  and  whether  it  is  a 
linear  or  non-linear  effect.  These  components  are  determined  in  a  manner  very  similar  to  that  used 
in  the  U.S.  Navy’s  Ship  Motion  Program  (SMP)  code  (Meyers,  et  al,  1981) .  However,  in  the  SMP 
code,  forces  are  calculated  in  the  frequency  domain,  assuming  certain  averaged  magnitudes  of  roll 
displacement  and  roll  velocity. 

Such  an  averaged  roll  damping  approach  is  not  satisfactory  for  time  domain  calculations  where  a 
primary  objective  is  the  accurate  calculation  of  the  extreme  response  events.  The  new  calculation 
method  uses  the  formulae  fi-om  the  references  in  Table  1,  but  uses  the  current  magnitude  of  roll 
displacement  and  roll  velocity  rather  than  an  averaged  value.  At  every  time  step,  the  time  history 
of  roll  displacement  and  roll  velocity  is  examined  for  a  peak  value,  positive  or  negative.  These 
peak  values  generate  parameters  for  the  viscous  forces  until  a  new  peak  is  found.  At  any  given  time 
step,  the  actual  forces  depend  on  these  parameters  and  the  instantaneous  value  of  roll  displacement 
and  roll  velocity.  This  approach  is  very  different  from  the  approach  used  in  SMP,  which  uses  an 
iterative  process  to  calculate  an  “equivalent”  or  “averaged”  roll  amplitude  for  viscous  damping. 
The  current  approach  is  a  more  direct  calculation  taking  advantage  of  the  fact  that  the  roll  angle  and 
velocity  are  known  at  all  time. 


Table  1.  Viscous  and  Lift  Effects 


Effect 

Reference 

Linearity 

Hull  Lift 

Himeno  (1981) 

Linear 

Skeg,  Bilge  Keel,  and  Foil  Lift 

Himeno  (1981) 

Linear 

Hull  Eddy-Making 

Tanaka  (1960)  and  Ideda  et  al.  (1978) 

Non-Linear 

Bilge  Keel  Eddy-Making 

Kato(1966) 

Non-Linear 

Skeg  and  Foil  Eddy-Making 

Tanaka  (1960),  Ideda  et  al.  (1978), 
and  Tamiya  (1972) 

Non-Linear 

Hull  Skin  Friction 

Kato(1966) 

Non-Linear 

2.2.2  Lift  Forces 


A  lift  force  is  calculated  for  every  foil,  including  the  hull  and  bilge  keels.  In  the  following 
discussion,  all  the  values  are  described  in  a  coordinate  system  fixed  to  the  ship.  Assume  we  have 
determined  the  ship  motion  and  the  fluid  flow  in  this  system: 

=  velocity  of  ship  at  local  origin 
=ship  rotation  rate  about  local  origin 

We  know  the  following  information  about  the  foil: 

5  =  ship  rotation  rate  about  local  origin 

C  =  average  chord  direction  fi-om  tail  to  nose 
P  foii  =  location  of  foil  center 

V  =  fluid  velocity  at  foil  center 

Then  we  can  determine  the  relative  velocity  at  the  foil  center,  as  follows: 

Vr^V-  (U  ship  +  'Rship  X  'Pfcii )  (1 ) 


We  normalize  these  three  vectors  as: 


•4 

5 


c  = 


VR 


To  find  the  velocity  for  the  lift  calculation,  we  determine  the  component  of  velocity  in  the  span 
direction,  and  subtract  it  from  the  total  relative  velocity. 

Ve=VR-S(VR-s)  (2) 


We  find  the  unit  flow  direction  as 
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(3) 


The  unit  lift  direction  is  perpendicular  to  this  effective  flow  and  the  span  direction, 

/  =  Ve  X  5  (4) 

The  effective  angle  of  attack  is  the  angle  between  the  flow  vector  and  the  chord  tail-to-nose  vector. 
This  is  a  rough  approximation  that  that  works  well  for  rectangular  shapes  and  is  sufficient  for  the 

purposes  of  this  formulation.  The  vectors  v^and  c  are  normalized,  so  the  cosine  of  the  angle  of 
attack  a  is  the  dot  product  of  the  two: 

ar  =  cos“  (ve  -5)  (5) 


The  lift  coefficient  is  determined  as  a  function  of  the  aspect  ratio  In  all  cases,  the  foils  are 
assumed  to  be  "groundboarded,"  with  the  effect  of  doubling  the  geometric  aspect  ratio.  The 
effective  aspect  ratio  is  then 


where:  S  =  effective  span  of  foil 

C  =  effective  chord  of  foil 


(6) 


This  is  the  result  of  Prandt’s  lifting  line  theory  and  is  sufficient  for  the  purposes  of  these 
calculations.  Then  the  lift-curve  slope  is: 


II 

2  " 

if  <  2.0 

Itt  , 

- 

1.7  ^ 

if  >  2.0 

(7) 

(8) 


The  foil  is  assumed  to  have  this  lift-curve  slope  up  to  a  limiting  angle  of  attack,  Above  this 
angle,  the  lift  goes  to  zero.  If  the  flow  is  at  a  small  angle  to  the  trailing  edge,  the  foil  is  assumed  to 
be  lifting  "backwards"  and  has  a  different  stall  angle, 

The  lift  coefficient  is: 


C,,  -a|F.|  (S  C) 

:  |a| 

C,,-{7t-aWe\  (S  C) 

:  |/r 

(9) 

0.0 

•  ^  stall  <  ^  ~  ^  stall 

0.0 

:  -;r 
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In  the  current  program,  the  stall  angles  are  as  follows: 


FoUs 

HuU 

Bilge  Keels 

24 

12 

12 

12 

12 

12 

To  properly  orient  the  lift  vector,  its  direction  relative  to  the  lifting  surface  must  be  the  same  as  the 
normal  component  of  the  incident  velocity.  To  determine  this,  we  take  the  sign  of  the  dot  product, 

Siift  =  signive  -  c(v  •  c)  •  /)  (1 0) 

where  S,^  =  the  direction  of  the  lift  influence 

The  lift  vector  is  then 

(11) 

where  p  =  density  of  water. 

In  the  current  version  of  LAMP,  the  velocity  on  the  appendages  does  not  include  the  wave  particle 
velocities,  but  can  include  a  velocity  increment  to  simulate  the  effects  of  the  propeller  wash. 
However,  this  velocity  increment  was  not  available  in  LAMP  when  these  calculations  were 
performed. 

2.2.3  Effect  of  Eddy-Making  Forces  on  Roll 

The  eddy-making  force  is  calculated  as  the  force  on  a  flat  plate  in  a  flow  normal  to  its  surface.  This 
is  the  time-domain  equivalent  of  the  semi-empirical  methods  used  in  SMP,  developed  by  Kato 
(1966). 

First  we  find  the  normal  flow,  Vn  ,  as 

V^=V-(cxs)  (12) 

The  Reynolds  number  based  on  the  span  is 

(13) 

V 

The  "coefficient  of  normal  pressure,"  C;^,  is  calculated  (after  Kato)  as 
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(14) 


C;^  =  1.98  exp- 1 1.0-^  :  <  0.07 

Cff  =  1.98  exp- 1  l.Of  :  f  <  0.07 

C^=1.18  :  0.07<|:<^ 

The  "normal  pressure  magnification,"  ,  is  calculated  (after  Kato)  as 


C^  =  1.95-0.251n(/?£)  +  0.20sin(^(ln(/?^)-2.19))  :  7?^  <1000  (15) 

Q  =  1.0  :  Re  ^  1000 

The  hull  form  correction,  C^,  and  the  factor  for  bilge  keels  built  up  fi’om  plates,  Q,  are  both  set  to 
1.0.  The  normal  or  eddy-making  force  is  then  calculated  as: 


(16) 


Any  drag  force  can  be  calculated  and  has  the  direction  of  the  incident  flow,  Ve.  In  the  current 
implementation,  the  flat-plate  skin  friction  drag  is  calculated  using  the  local  Reynolds  number. 


2.2.4  Effect  of  Hull  Skin  Friction  on  Roll 

The  following  derivation  shows  how  the  semi-empirical  frequency  domain  techniques  are  applied 
to  the  time-domain  problem.  Kato  (1958)  established  a  semi-empirical  technique  for  estimating  this 
component.  His  formulae  expressed  the  damping  as  a  roll  decrement  ratio  for  zero  forward  speed. 
Further  work  by  Tamiya  (1972)  added  the  effect  of  speed,  and  Himeno  summarized  this  work  in 
1981,  expressing  the  damping  coefficient  as  follows: 

+  (17) 

CQL  2 

where  644^. 

U 
L 

CO 

P 
S 


^4(0 

^4(0 


=  the  roll  damping  due  to  skin  fnction  with  forward  speed 
=  ship  speed 
=  ship  length 
=  roll  frequency 
=  water  density 
=  surface  area 
=  "equivalent"  radius 
=  fnctional  coefficient 

=  roll  amplitude 
=  roll  velocity  amplitude 
=  time-dependent  roll  angle 
=  time-dependent  roll  velocity 
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We  write  the  frequency-domain  moment  amplitude  as  the  product  of  the  damping  coefficient  and 
the  roll  velocity  magnitude: 

M^y=(\+A.\^){]-pSr^Cfr]^^)  (18) 

COL  2 

In  the  quasi-steady  approach,  we  vary  both  the  damping  coefficient  and  the  roll  velocity  with  tme. 
To  use  the  frequency-domain  data,  we  need  to  choose  an  appropriate  value  for  the  magnitudes 
and  f]^ .  At  every  time  step,  we  check  the  value  of  Ti^{t)  and  to  see  if  they  have  passed 
through  zero.  If  they  have,  we  take  the  previous  maximum  as  the  new  value  of  the  roll  magnitude. 
At  every  time  step,  we  calculate  the  moment  as 

(19) 

COL  2 

This  is  the  procedure  followed  at  every  time  step  i  to  track  the  value  of  the  roll  magnitudes: 

1 .  If  the  sign  of  does  not  equal  the  sign  of  774, find  7/4  as  the  maximum  since  the  last  zero 

crossing  _ 

2.  If  the  sign  of  ^4,.  does  not  equal  the  sign  of  7/4, find  ^4  as  the  maximum  since  the  last  zero 

crossing 

3.  Calculate 

2.3  Seaway  Description 

2.3.1  Overview 

The  seaway  definition  used  in  the  LAMP  system  allows  any  discrete  set  of  waves  to  be 
superimposed.  The  program  LAMP  seaway  module  allows  the  user  to  generate  one  of  these  sets 
based  on  parameters  using  standard  spectral  seaway  definitions  found  from  empirical  data.  The  user 
defines  the  power  spectrum  of  the  seaway  with  either  the  Brettschneider  Two-Parameter  Spectrum, 
which  uses  the  significant  wave  height  and  the  modal  period,  or  the  Pierson-Moskowitz  one- 
parameter  spectrum,  which  uses  only  the  significant  wave  height. 

2.3.2  Brettschneider  Two-Parameter  Spectrum 

The  user  supplies  the  modal  period  and  significant  wave  height  of  the  seaway. 

The  Brettschneider  Power  Spectrum  is  evaluated  as  follows: 

S(Q))  =  486.00* (20) 
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where 


Q)  =  wave  frequency 

2^  =  modal  period  (wave  period  at  which  S(^(0)  is  maximum) 

=  significant  wave  height 

3 

S{co)  =  amplitude  spectrum  for  the  seaway. 


2.3.3  Pierson-Moskowitz  One-Parameter  Spectrum 


If  a  Pierson-Moskowitz  spectrum  is  chosen,  only  the  significant  wave  height  is  input,  and  the 
modal  period  is  calculated  as  below.  Note  that  the  Pierson-Moskowitz  formula  is  generally  based 
on  wind  speed,  but  that  there  assumes  the  following  direct  relationship  between  wind  speed  and 
significant  wave  height. 

fj  -El.  / 16.034*0.0081'  „!) 

1  "  g  V  4*0.74 

where  U  is  the  wind  speed  10  meters  above  the  free  surface. 


Then,  if  is  either  input  or  calculated  based  on  U ,  the  modal  period  2^  can  be  calculated  as 

3 

follows: 


°  g/ci  a 

where  g  =  gravitational  acceleration 

k  =  4.00043 

a  =  0.0081,  Phillips  Constant. 


(22) 


2.3.4  Discretization  and  Simulation  of  a  Time-Domain  Wave 


In  order  to  discretize  the  spectrum,  the  program  has  the  user  input  the  number  of  waves  to  use  and 
chooses  the  number  of  waves  to  define  the  seaway.  The  frequencies  are  determined  by  a  geometric 
spreading  away  from  the  modal  frequency  in  both  directions.  Half  of  the  specified  numbers  are 
used  for  frequencies  below  the  modal  frequency,  and  half  above.  We  find  the  lowest  frequency  as 
that  which  has  spectral  energy  not  greater  than  0. 1  percent  of  the  maximum  spectral  energy, 


1000. 


(23) 


Then  we  determine  geometric  spacing  for  N,  points  from  to  Wq  as  follows: 

<^=1.25 


0),  =0} 

*^min 

A  ^0 

ACO:  =  - ^ 


■0)„ 


^  -  0.5(1 +  (J"-') 
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(24) 

(25) 

(26) 


The  "central"  frequencies  are  then 


ffl,.  =  +  Oi*  ,  for  /  =  2  to  iV, 


(27) 


We  distribute  N2  =  N-N^  frequencies  above  (Oq  in  the  same  way,  this  time  choosing  a  value  of  S 
that  yields: 


1000. 


(28) 


Then  the  amplitudes  of  the  N  waves  representing  the  spectrum  are  found  as 

^,=V25(ty,)Aty,  (29) 

The  phases  are  found  using  a  uniform  random  distribution  from  ;r  to  -;r. 


2.4  Simple  PID  Rudder  Control  Algorithm  in  LAMP 


When  a  ship  is  operating  in  oblique  seas,  it  is  necessary  to  have  rudder  control  to  keep  the  ship  on 
course.  This  is  routine  in  model  tank  experiments.  In  the  current  version  of  the  LAMP  code,  a 
simple  PID  (Proportional,  Integral,  and  Derivative)  control  algorithm  in  added  for  course-keeping 
during  ship  motion  simulation.  Based  on  the  current  and  desired  heading  angles  of  the  ship,  the 
following  rudder  command  is  given: 


where 


=  rudder  command 

y/  =  ship  heading  at  the  current  time  step 

y/^  =  ship  heading  at  the  previous  time  step 

y/^  =  desired  ship  heading 

=  gain  coefficient  for  the  proportional  term 
=  gain  coefficient  for  the  derivative  term 
At  =  time  step  size. 


(30) 


Only  proportional  and  derivative  terms  are  used.  The  integral  term  that  is  usually  used  to  correct 
certain  biases  is  not  included  in  this  formulation.  In  the  current  control  algorithm,  G„  is  set  to  be 
0.9  and  G^  is  set  to  be  20.0. 

For  rudder  dynamics,  it  is  assumed  that  the  rudder  servo  is  a  linear  first  order  lag  with  given  slew 
rate  R„,  and  proportional  band  Pg.  Thus, 

t5+6  =  5,  (31) 
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where  t=  P^l  In  the  current  algorithm,  Pg  is  set  to  be  5  degrees  and  is  set  to  be  5.5 
degrees/second.  If  the  rudder  command  6^  is  constant  in  the  time  interval  At,  the  exact  solution  of 
equation  is 

S=5,  +  {\-e-"){8^-5,)  (32) 


where  5q  is  the  rudder  angle  at  the  beginning  of  the  time  step  and  the  second  term  is  the  increment 
of  the  rudder  angle.  Note  that  this  increment  is  limited  by  the  rudder  slew  rate  i.e. 


At 


<R^ 


(33) 


The  rudder  angle  is  limited  by  the  maximum  rudder  angle  ,  i.e. 


(34) 


In  the  current  control  algorithm,  is  set  to  be  35  degrees. 


2.5  A  Multi-Level  System 

A  complete  computational  capability  for  the  assessment  of  ship  motions  and  wave  loads  must  be 
based  on  a  multi-level  approach.  Such  a  system  integrates  mediods  that  are  based  not  just  on  one 
single  code  or  one  single  level  of  sophistication,  but  rather  on  a  system  of  codes  with  different 
levels  of  sophistication.  As  a  general  rule,  the  physics  underlying  the  ship/wave  interactions  is  best 
understood  using  comparisons  generated  by  incremental  increases  in  complexity  -  a  procedure  that 
also  moderates  computer  usage.  Analysis  tools  at  the  lower  levels  may  employ  several 
approximations  to  attain  a  short  enough  turnaround  time  for  use  in  early  stages  of  the  evaluation 
process.  Examination  of  results  obtained  by  the  lower  level  code  guides  the  engineer  in  choosing 
areas  where  more  accurate  theories  must  be  used.  In  other  words,  the  lower  level  codes  should  be 
used  as  a  filtering  mechanism  for  the  selection  of  more  accurate  but  more  complicated  and 
computationally  intensive  codes. 

A  multi-level  system  can  also  effectively  tie  the  probabilistic  and  deterministic  approaches  together 
providing  the  missing  ingredient  of  probabilistic  prediction.  Statistical  data  of  ship  motion  in  given 
random  seas  can  be  obtained  by  using  lower  level  evaluation  codes  to  efficiently  compute  the  ships 
responses  to  a  very  wide  range  of  deterministic  excitations.  The  severe  ship  responses  can  be 
selected  from  these,  to  be  examined  with  the  higher  level  nonlinear  simulations.  Conversely, 
nonlinear  dynamic  simulations  of  ships  in  episodic  wave  events  can  be  used  to  understand  the 
actual  physical  mechanisms  underl3dng  the  ship  responses  to  these  events,  such  as  capsizing,  and  to 
identify  dominant  factors  of  vessel  stability,  which  can  be  used  in  the  statistical  screening  process 
using  the  lower  level  codes. 

Recognizing  these  needs,  the  LAMP  System  is  being  developed  as  a  multi-level  code  system 
consisting  of  a  total  of  three  computational  methods  of  different  levels  of  sophistication. 

•  LAMP-4:  The  large-amplitude  3-D  nonlinear  method 

•  LAMP-2:  The  approximate  large-amplitude  3-D  nonlinear  method 

•  LAMP-1 :  The  linearized  3-D  time-domain  method 
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The  LAMP-4  method  is  the  complete  large-amplitude  method  where  the  3-D  potential  is  computed 
with  the  linearized  free-surface  condition  satisfied  on  the  surface  of  the  incident  wave.  Both  the 
hydrodynamic  and  hydrostatic  pressures  are  computed  over  the  instantaneous  hull  surface  below 
the  incident  wave  surface.  Large  computer  resources  are  required  for  this  method.  In  the  LAMP-2 
method,  the  linear  3-D  approach  is  used  to  compute  the  hydrodynamic  part  of  the  pressure  forces, 
while  the  hydrostatic  restoring  and  Froude-Krylov  forces  are  calculated  with  the  same  accuracy  as 
in  LAMP-4.  The  reason  for  developing  this  simplified  method  is  that  it  drastically  reduces  the 
requirements  for  computer  resources.  The  LAMP-1  method  is  the  linearized  version  of  the  LAMP- 
4  method.  This  3-D  time-domain  method  includes  a  routine  for  automatic  generation  of  the 
fi"equency  domain  results. 

Table  2  shows  how  the  hydrostatic  restoring  and  Froude-Krylov  forces  and  the  hydrodynamic 
(added  mass,  damping  and  diffraction)  forces  are  calculated  for  the  four  different  LAMP  methods. 
The  hardware  requirements  for  the  four  methods  are  also  shown  in  Table  2.  Note  that  the  two 
nonlinear  methods,  LAMP-2  and  LAMP-4  are  based  on  the  approach  that  both  the  motions  and  the 
waves  may  have  large  amplitudes.  For  all  of  these  three  nonlinear  methods,  the  restoring  and 
Froude-Krylov  forces  are  calculated  exactly  over  the  instantaneous  wetted  surface  below  the 
incoming  wave  surface. 


Table  2.  Computation  Methods  and  Hardware  Requirements  for  the  LAMP  Code 

(Z=0  and  F(t)  are  Still  Water  Surface  and  Incident  Wave  Surface,  Respectively) 


Method 

Hydrodynamic,  Restoring,  and  Froude-Krylov  Forces 

Hardware 

LAMP-4 

Free  Surface  Boundary  Conditions  on  F(t) 

3-D  Large-Amplitude  Hydrodynamics 

Nonlinear  Restoring  and  Froude-Krylov  Forces 

Fast  Workstation 

LAMP-2 

Free  Surface  Boundary  Conditions  on  F(t) 

3-D  Linear  Hydrodynamics 

Nonlinear  Restoring  and  Froude-Krylov  Forces 

Workstation 

LAMP-1 

Free  Surface  Boundary  Conditions  on  Z=0 

3-D  Linear  Hydrodynamics 

Linear  Restoring  and  Froude-Krylov  Forces 

Workstation 
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3  Fishing  Vessel  Model  and  Static  Stability 

3.1  Model  for  LAMP 


Figure  3  shows  a  body  plan  view  of  the  fishing  vessel  Italian  Gold ,  as  it  was  designed.  As  noted  in 
the  Coast  Guard  Report  (USCG,  1994),  these  plans  do  not  correspond  exactly  to  the  vessel  as  built. 
These  design  drawings  are  the  best  available  model  of  the  hull,  and  were  used  by  the  Coast  Guard 
in  their  analysis. 


Figure  3.  Body  Plan  of  the  Fishing  Vessel 


The  LAMP  program  requires  a  full  surface  definition,  represented  by  an  array  of  quadrilateral 
panels.  This  geometry  was  created  by  measuring  the  ship’s  molded  keel  line,  chine  line,  and  main 
and  whaleback  deck  lines.  A  reasonable  approximation  to  the  developable  surface  was  obtained  by 
distributing  points  along  each  of  these  lines  in  proportion  to  arc  length.  Figure  4  shows  this 
panelization  in  a  body  plan  view.  The  left  side  shows  the  original  panelization.  In  order  to  reduce 
the  number  of  surfaces  for  the  LAMP  calculation,  the  chine  was  "softened,"  so  that  the  lower  hull 
and  the  upper  hull  could  be  joined  as  one  surface  without  a  hard  line.  This  geometry  is  shown  on 
the  right  side  of  Figure  4.  A  profile  view  of  the  "softened"  geometry  is  shown  in  Figure  5. 


16 


X(Ft.) 


Figure  5.  Profile  View  of  Softened  Panelization 


3.2  Hydrostatic  Stability 

An  evaluation  of  the  undamaged  hydrostatic  stability  of  the  vessel  was  performed  using  LAMP 
hydrostatics  and  data  from  the  Coast  Guard  Report  on  the  casualty  (USCG,  1994).  This  report 
evaluated  the  vessel  in  its  designed  configuration,  not  "as-built".  The  design  drawings  (Gilbert 
Associates,  1979)  were  used  for  the  USCG  geometry  as  well  as  the  LAMP  geometry  for  an 
assumed  loading  condition  as  follows: 


Loading  Condition 

English  Units 

Metric  Units 

Displacement 

190.73  LT 

193.79  MT 

VCG  above  Baseline 

8.87  ft 

2.70  m 

LCG  Relative  to  Amidships 

2.77  ft  Aft 

0.88  m 
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The  LAMP  hydrodynamic  model  did  not  include  the  skeg,  the  rabbet,  the  propeller,  the  rudder  and 
other  appendages.  These  appendages  were  included  as  components  generating  viscous,  but  not 
hydrodynamic  forces.  The  hydrostatics  generated  by  LAMP  reflect  these  differences.  The  intent  of 
the  LAMP  hydrostatic  calculations  was  to  establish  the  static  stability  of  the  vessel  as  seen  by 
LAMP  for  our  calculation.  When  applying  the  displacement  of  190.73  Long  Tons  (193.79  MT),  the 
LAMP  hydrostatics  gave  a  baseline  draft  amidships  of  8.422  ft  (2.567  m),  deeper  than  the  USCG 
draft.  The  USCG  report  did  not  provide  the  draft  at  zero  heel  and  trim  for  the  undamaged  condition, 
but  the  baseline  draft  for  heel  angle  of  3.89  degrees  was  given  as  7.827  ft  (2.386  m).  We  chose  to 
match  the  displacement  rather  than  the  draft. 

LAMP  was  used  to  generate  the  vessel’s  righting  arm  or  GZ  curve  versus  heel  angle.  The  vessel  is 
given  a  fixed  heel  angle  and  lowered  into  the  water.  When  the  vertical  buoyancy  force  matches  the 
displacement,  the  transverse  moment  about  the  center  of  gravity  is  measured.  This  quantity  is 
divided  by  the  displacement  to  give  the  righting  arm. 

In  Figure  6,  the  LAMP  GZ  curve  is  plotted  along  with  the  wind-heeling  arm  from  the  USCG 
calculations.  The  wind  force  is  based  on  a  53.4  knot  steady  wind  abeam.  It  takes  into  accoimt  the 
components  of  the  above-water  structure.  The  first  intersection  of  these  curves  shows  the  static  heel 
angle  that  would  be  produced  by  this  wind  abeam,  approximately  6  degrees.  The  second 
intersection  shows  the  maximum  angle  due  to  wave  action  that  can  be  sustained  before  capsize  in 
this  wind,  approximately  46  degrees.  This  analysis  assumes  a  calm  free  surface,  while  the  actual 
hydrostatic  stability  in  a  seaway  will  depend  on  the  actual  wave  elevation. 


CO 


Figure  6.  Static  Stability  of  Fishing  Vessel  with  USCG  Wind  Heel 
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4  Nonlinear  Motion  Simulation  of  Fishing  Vessel  in 
Stern  Quartering  Seas 

To  demonstrate  the  application  of  the  improved  LAMP  code  for  assessing  the  dynamic  stability  of 
a  fishing  vessel  in  waves,  a  series  of  computations  were  performed  using  LAMP-2  for  the  Italian 
Gold  in  regular  and  irregular  stem  quartering  waves.  Only  intact  stability  is  studied.  The  center  of 
gravity,  CG,  of  the  ship  is  located  amidship  and  is  about  0.27  ft  (82.29  cm)  above  the  design 
waterline  in  the  current  study.  The  actual  location  of  the  center  of  gravity  depends  on  the  weight 
distribution  of  the  ship  and  is  a  very  important  factor  in  the  roll  stability.  The  CG  used  in  the 
current  calculations  is  a  reasonable  one  for  such  a  fishing  vessel  in  full  load  condition. 

In  principle,  the  LAMP  code  is  applicable  to  general  large-amplitude  motion  simulation,  including 
capsizing,  of  ships  in  a  wide  range  of  sea  conditions.  For  the  current  study,  uni-directional  linear 
incident  waves  are  selected.  The  wave  conditions  for  the  runs  with  regular  incident  waves  are  listed 
in  Table  3  and  for  the  runs  with  irregular  incident  waves  represented  by  Pierson-Moskowitz 
spectrum  are  listed  in  Table  4.  Note  that  n)is  the  wave  frequency,  A.  is  the  wavelength,  h  is  the 
wave  height,  A1/3  is  the  significant  wave  height,  and  L  is  the  ship  length  (68  ft,  20.7  m).  T^e  wave 
amplitude  and  wave  frequency  are  varied  to  see  the  effects  of  these  parameters  on  capsizing.  All 
waves  selected  are  bounded  by  the  "steepest  wave"  limit. 


Table  3.  Regular  Wave  Cases  of  LAMP  Runs 


Case  No. 

(rad/sec) 

A.  (ft) 

H(ft) 

X  (m) 

h(m) 

h/L 

hIX 

R-1 

1.7000 

70 

6 

21.3 

1.8 

0.088 

0.086 

R-2 

1.4224 

100 

6 

30.5 

1.8 

0.088 

0.060 

R-3 

1.2021 

140 

6 

42.7 

1.8 

0.088 

0.043 

R-4 

0.8500 

280 

6 

85.3 

1.8 

0.088 

0.021 

Table  4.  Random  Wave  Cases  of  LAMP  Runs 
One  Parameter  Pierson-Moskowitz  Spectrum 


Case  No. 

/Jl/3(ft) 

hm  (m) 

IR-1 

10 

3.1 

lR-2 

13 

4.0 

IR-3 

15 

4.6 

As  for  the  ship  operating  condition,  the  forward  speed  is  set  to  be  10  knots;  the  initial  roll,  pitch, 
and  yaw  angles  are  zero;  and  the  ship  is  at  the  bottom  of  the  incident  wave.  The  incident  wave  is 
coming  from  the  starboard  direction  of  the  stem.  The  angle  between  the  direction  of  the  wave 
propagation  and  the  direction  of  the  ship  forward  motion  is  45  degrees.  In  all  simulation  mns,  a 
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PID  control  is  used  for  course-keeping  purpose.  Six-degree-of-freedom  coupled  equations  of 
motions  are  solved  at  each  time  step  for  motion  simulation. 

A  sequence  of  snap  shots  for  case  R-2  (X,=100  ft,  30.5  m)  is  given  in  Figure  7.  These  snap  shots  are 
approximately  one  second  apart  during  a  10-second  simulation  run.  The  ship  started  on  course  and 
the  bow  was  pushed  toward  the  port  side  by  the  stem  quartering  wave  coming  from  the  starboard. 
The  autopilot  system  tried  to  maintain  the  original  course  by  turning  mdder  toward  the  port  side. 

As  can  be  seen  in  the  picture,  the  ship  bow  was  turning  back  and  at  the  same  time  rolling  toward 
the  port  side.  The  second  wave  crest  hit  the  ship  from  the  starboard  at  the  time  it  was  rolling  toward 
the  port  side.  As  a  result,  the  ship  capsized  due  to  dynamic  effect.  This  is  a  typical  broaching 
phenomenon.  The  behavior  of  the  ship  is  closely  related  to  the  wave  characteristics. 

Plots  of  the  motion  histories  in  six  directions  and  snap  shots  for  the  plotted  case  are  given  from 
Figure  8  to  Figure  1 1.  In  these  figures,  the  dimension  for  linear  displacement  is  ft  and  the 
dimension  for  the  rotation  is  degrees.  Notice  that  when  the  wavelength  is  about  1  to  2  ship  lengths, 
the  ship  will  capsize  even  though  the  wave  height  is  only  6  ft  (1.8  m).  When  the  wavelength  is  four 
times  the  ship  length  (case  R-4,  Figure  1 1),  the  ship  is  just  riding  on  the  wave  and  does  not  capsize. 

Numerical  simulations  have  also  been  done  for  Italian  Gold  in  fiilly  developed  random  seas.  The 
motion  histories  for  three  cases  (IR-1  with  hi/3=\0  ft  (3. 1  m),  IR-2  with  Ai/i=13  ft  (4.0  m),  and  IR-3 
with  hi/f=\S  ft  (4.6  m)  are  shown  in  Figure  12  to  Figure  14.  It  can  be  seen  clearly  that  the  ship 
capsized  while  the  significant  wave  height  reached  13  ft  (4.0  m).  It  is  interesting  to  see  that  the  ship 
can  survive  in  10  ft  (3.1  m)  random  sea  but  will  capsize  in  6  ft  (1.8  m)  regular  sea.  It  also  indicates 
that  the  most  dangerous  condition  in  a  random  sea  is  when  several  waves  group  together  like  a 
regular  sea  condition. 
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Figure  7.  Case  R-2:  Snap  Shots  of  Motion  Animation  of  Fishing  Vessel  Italian  Gold  in 
Linear  Regular  Stem  Quartering  Waves  with  ?i=100  ft  (30.5  m) 
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Figure  8.  Case  R-1:  Time  History  of  Motions  in  Six  Directions  for  Fishing  Vessel 
Italian  Gold  in  Linear  Regular  Stern  Quartering  Waves  with  X,=70  ft  (21.3  m) 
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Figure  10.  Case  R-3:  Time  History  of  Motions  in  Six  Directions  for  Fishing  Vessel 
Italian  Gold  in  Linear  Regular  Stern  Quartering  Waves  with  A,=140  ft  (42.7  m) 
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Figure  11.  Case  R-4:  Time  History  of  Motions  in  Six  Directions  for  Fishing  Vessel 
Italian  Gold  in  Linear  Regular  Stern  Quartering  Waves  with  X=280  ft  (85.3  m) 
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5  Summary 


A  new  numerical  simulation  method,  LAMP,  has  been  developed  for  studying  the  extreme  motion 
including  capsizing  of  ships  in  oblique  seas.  A  sample  study  of  the  static  and  dynamic  stability,  of  a 
typical  fishing  vessel  Italian  Gold  is  presented  in  this  report.  LAMP  motioii  simulations  of  the  ship 
in  various  incident  waves  were  performed.  It  was  found  that  under  the  specified  loading  condition, 
the  ship  will  capsize  in  stem  quartering  sea  condition  while  the  wavelength  is  about  1  to  2  times  the 
ship  length  and  wave  height  is  6  ft  (1 .8  m). 

On  the  other  hand,  the  ship  is  able  to  survive  10  ft  (3.0  m)  significant  wave  height  in  a  fully 
developed  random  sea.  However,  the  condition  reported  at  the  time  while  Italian  Gold  capsized 
was  55  knots  NE  wind.  This  corresponds  to  sea  state  7  with  significant  wave  height  of  19  ft  (5.8 
m).  As  shown  in  the  numerical  simulation,  the  ship  would  not  be  able  to  survive  at  the  given 
loading  condition. 


Only  intact  stability  was  considered  in  this  study.  From  the  numerical  results,  it  is  found  that  water 
on  deck  may  be  important  for  ships  in  several  different  wave  conditions.  Loads  due  to  water  on 
deck  or  possible  compartment  flooding  should  be  taken  into  consideration  in  the  future  study.  Wind 
force  may  be  another  important  factor.  Wind  force  will  create  steady  heeling  moment,  which  will 
further  reduce  the  stability.  Both  of  these  factors  are  unfavorable  to  the  stability  of  the  ship. 

The  current  LAMP  simulation  is  full  six-degree-of-ffeedom  with  speed  condition.  Although 
viscous  and  lifting  effects  are  included  (skin  fnction,  bilge  keel  and  rudder).  Several  other 
important  factors,  such  as  the  propeller  thrust,  the  effect  of  the  propeller  slipstream  on  the  mdder, 
the  effect  of  boundary  layer  separation  at  the  stem  of  the  vessel  on  the  maneuvering  characteristics 
of  the  vessel,  wind  effects,  and  nonlinear  wave  effects,  were  not  modeled  in  this  simulation  and 
should  be  included  in  any  future  studies. 
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